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Abstract 

We present the results of numerical experiments in connection with the Riesz and Hardy- 
Littlewood criteria for the truth of the Riemann Hypothesis (RH). The coefficients Ck of the 
Pochammer's expansion for the reciprocal of the Riemann Zeta function, as well as the "critical 
functions" Ckk a (where a is some constant), are analyzed at relatively large values of k. It appears 
an oscillatory behaviour (Riesz-Hardy-Littlewood wave). The amplitudes and the wavelength of the 
wave are compared with an analytical treatment concerning the wave in the asymptotic region. The 
agreement is satisfactory. We then find numerically that in the large (5 limit too, the amplitudes of 
the waves appear to be bounded. For a special case the numerical experiments are performed up to 
larger values of k, i.e k = 10 9 and more. The analysis suggests that RH may barely be true and an 
absolute bound for the amplitudes of the waves in all cases should be given by | ^i +t ^ ~M> with e 
arbitrarily small positive, i.e. equal to 1.68. . . . 
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I. INTRODUCTION 



Following recent works concerning the study of some well known functions appearing in 
the original criteria of Riesz, Hardy and Littlewood for the possible truth of the Riemann 
Hypothesis (RH), there is new interest in the direction of numerical experiments, where the 
calculations use the ideas of some recent works on the subject. These concern the expansion 
of the reciprocal of the Riemann Zeta function in terms of the so called Pochammer's poly- 
nomials Ph, whose coefficients play a central role also in the asymptotic region of very 
large k Q, Q, [], U]. For new zero free regions of the Zeta function, in the context of a rigorous 
treatment with the Miintz formula, the reader may consult a recent work by Albeverio and 
Cebulla fl. 

Here we are concerned with the discrete version of the Riesz criterion which has also 
been studied numerically: the first numerical experiments for values of k up to lOO'OOO 
have been announced and reported for the Riesz case in j^, [jj. It has been found that 
the function has a oscillatory behaviour in a region of relatively high fc's, in agreement 
with an asymptotic formula given by Baez-Duarte |2|. The agreement appears satisfactory 
even if only the contribution of the first non trivial zero of the Zeta function located at 
s = \ + % 14. 134725 in the complex plane has been used. 

In a previous work |7] a two parameter family (parameter a and [3) of Pochammer's 
polynomials was introduced. This allowed the starting investigation of at low values of k, 
but in various cases and in the so called "strong coupling" regime (high (3). After the initial 
study at low k, our computations using the formula containing the Mobius function were 
easily extended to larger and larger k (up to a billion) in the strong coupling limit, with 
the appearance of macro-oscillations in Cfc extending to larger k. This is a symptom that 
using such a limit the RH may eventually barely be trueQ. In this work we continue the 
numerical experiments also using our Poisson formula established in Q which is well suited 
for numerical purposes. 

After the formulation of the model in Section 2, we then compute in Section 3 the 
amplitudes of what we call the Riesz-Hardy-Littlewood wave, which contains arbitrarily 
scales in few of the two parameters a and (3 at our disposal. Using the Baez-Duarte formula, 
we then present our results for different models up to values of k equal to one billion and 
observe oscillations in all cases (Section 4). The agreement with the asymptotic formula of 
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Baez-Duarte is satisfactory. Then, in Section 5, we concentrate the study in more details 
by considering a special new model already proposed in Q where a — | and (3 is increasing 
starting with the value equal to 4. The results show in a concrete way the "transition" from 
the low coupling to the "strong coupling regime": at low values of (3 ((3 = 4) we obtain up to 
7 oscillation with values of k extending up to a billion. These start to deform continuously 
with increasing values of (3 approaching the infinite (3 limit. In such a regime, the wave 
is absorbed in a macroscopic region with an amplitude whose strength should be finite as 
already noted in 

In the context of validity of our numerical results, our analysis gives further indication 
that the RH may barely be true as indicated by our two parameter models in the week as 
well as in the "strong coupling regime" (Section 5). Moreover, the possibility that in an ideal 
numerical experiment (using an arbitrarily large but finite maximum value of n, say N in 
the formula with the Mobius function) the amplitude of the waves at finite (3 values should 
decrease, is commented in Appendix. 



II. THE MODEL 

Following recent treatments [l|, Q- Q], a possible expansion of the reciprocal of the Zeta 
function i.e. ((s) 1 in terms of the so called Pochammer's polynomials P^, with two param- 
eters a (a > 1) and (3 is this one: 



77~T = c fc( a ' P) p k(. s > oi, (3) 
C(s) z - 



k=0 



where 



— + 1 



P k (s,a,P) = l[ 1 

r=l \ 

n=l x 

and Pfc(0, a, (3) = 1. 

In (3) the Mobius function of argument n is given by: 



1, if re = 1 

(— l) fc , if re is a product of k distinct primes 
0, if re contains a square 



(1) 

(2) 
(3) 



If s = p + it is a complex variable and argument of the Riemann Zeta function one has 
for JR(s) = p > 1: 



_l_ = ^Mn) 

C(s) ^ n s 1 J 



n=l 

Another explicit formula for the Ck(a,0) is obtained from (3) using the binomial coeffi- 
cients and reads: 



As /3 is increasing, one may a 



the formula recently obtained 
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3=0 

so use (especially) in the context of numerical experiments, 
and given by: 

n=l 



In such an approximation we have that 

c fc (a, /?) = V c p (a, /?)-re- fc (7) 

which shows the emergence of a Poisson like distribution for the coefficients Ck(cc,(3). This 
should be a very satisfactory approximation in the limit of relative large values of j3 pf. We 
recall that an important inequality due to Baez-Duarte concerning the Pochammer's 
polynomials of complex argument z is given by: 

\Pk{z)\<Ck-^ (8) 

The above inequality applied to our two parameter family of Pochammer's polynomials 
with complex argument z = + 1 gives: 

\P k {s,a,i3)\<Ck-^ +1 (9) 

So that in (1) will be different from zero and thus the RH will be true for 9t(s) = p > \ 
if the Ck decay, at large k, as (see Q|): 

\c k \ < Ak~^ (10) 
We will also consider the "critical function" 

iP(k;a,{J,p):=c k k^ (11) 



TABLE I: The expected decay of c& for different values of a and ft 



a 


ft p decay of c& 


Note 


2 


o 1 
2 2 


k 4 


The case of Riesz 


1 


2 2 


fc-z 

K 4 


The case of Hardy-Littlewood 


2 


1 


1 


Same decay as the Hardy-Littlewood case but numerically more convenient 


7 
2 


4 i 
* 2 




Ddlllc LI.cL.d_y do lllc XVlcoZ L,dbc ; llllcIlblVc L,dlL.U.ldllU.lJ.o die Elvcll UclUW 


3 


3 i 
2 




o 3 

If the Zeta function has no zero for p > ^ then Cfc(3, 3) should decay at least as k * 


4 


4 i 

H 2 




Since from the Prime number theorem there is no zero for p > 1 
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the Cfc(4, 4) decays at least as fc~4 


2 


4 i 

* 2 




Another interesting case for calculations 



which from (10) is expected to be bounded by a constant A. 

We now recall two original cases given in pionnering works by Riesz |S| and by Hardy 
Littlewood Setting p = ~ in (10), for a = ft = 2 (Riesz case) we have that |c&| < Ak~% 
and for a = 1, ft = 2 (Hardy-Littlewood case) |c&| < Ak~i. Other interesting cases for which 
we will carry out intensive numerical experiments to be presented below are summarized in 
the Table I. 

A limiting delicate case analyzed in [3] is the one where a = \ + 5 and ft grows to 
infinity. Here of course we do not have absolute convergence to 1 (ck may nevertheless 
be analyzed) and from (10) we have that the c& should be smaller then a constant for all k. 
This is what we verified with numerical experiments (not presented here) with values of k 
up to a billion. The value of the constant has been proposed in our previous work f7j and 
the conjecture was that |c&| < I77T7 — 1| = 1.68477. However the situation is delicate (a < 1) 

n (5) - 

since Littlewood has shown that, assuming RH is true, ^ -^tt- is convergent for all e 

n=i n7 

strictly greater than zero. 

The general situation is that the "critical function" cyk~ should be bounded by a con- 
stant in absolute value as k — ► oo. In fact the function starts at zero for k = 0, reaches a 
minimum, then starts to increase and then begins to oscillate with a "constant " amplitude 
as k — > oo as we will see in the experiments. In a previous work pj we have analyzed c& in 
various cases but only for moderately values of k, i.e for k not exceeding 1000, with exception 
of some cases at large values of ft, where k reached the value of a half billion. Ck was found 
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to have only negative values in the range considered and increasing with k. Presently we 
know of recent numerical experiments ja| in the Riesz case carried out by Maslanka (k up 
to lOO'OOO) and Wolf (k up to 200'000). These calculations show that Ck become of oscil- 
latory type, thus assuming positive and negative values with an amplitude which appears 
constant in the range considered. In fact two or three oscillations with a wavelength related 
in first approximation to the first zero of the Riemann Zeta function may be seen. Here it 
should be remarked that this situation for the Riesz case is not in contraddiction with our 
strong coupling limit (0 large) cited above (see discussion below for the case a — ~ and (3 
increasing). 

In few of these new finding, we want first analyze (in an analytical context) such a 
behaviour and we call this general phenomena the Riesz-Hardy-Littlewood wave. This will 
be analyzed using an interesting result of Baez-Duarte, i.e. an expression giving c& for 
k — > oo. 



III. THE RIESZ-HARDY-LITTLEWOOD WAVE 

For the Riesz case, in connection with the Mellin inversion formula, the Riesz function is 
given (see Q] and (ll|) explicitly by: 

Using the calculus of residues F(x) is obtained by an integration and is given by: 

a+ioo 

, ^ i f r(l - s)x s , , , 



a—ioo 



where ~ < a < 1. 



Now, recently Baez-Duarte [2|], with an ingenious method found in particular an expres- 
sion for the reciprocal of the Pochammer polynomial given by: 



p k(s) \jjs-j 



where uniformely on compact subsets one has: 



6 



and he was able to obtain an explicit formula connecting Ck and the set of all trivial and non 
trivial zeros (z denote the complex Zeta zeros, z — \ + it) under the assumption of simple 
zeros. For the Riesz case the expression is given by: 

-^- = gcwW (16) 

for sufficiently large k. It should be said that formula (16) of Baez-Duarte is very nice and 
may be used to control our numerical computations at large k to be presented below. Ap- 
parently (16), with some precautions, may be extended to the general case with parameters 
a, (3 and should read: 

-^-'^ cwniW') (17) 

To obtain an asymptotic value for the amplitude of the Riesz-Hardy-Littlewood wave, we 
use (15) in (17): 

7 ^7 3^1 

3(z) ^ ^ ' 

where t = 3(z). In the limit of large k, one may neglect the contribution of the trivial zeros 
P|. For the "critical function" we then have the following expression: 

-i ir,* f r(-to _ 1 

k—c k = — 3 Y, c(z) =■ ^ ^ 2 } (19) 

for large k. To prepare the comparison of (19) with the numerical results we write explicitly 
(19) for the various cases we will treat. In order to obtain an estimate for the amplitude 
of the wave in the long wavelength limit (k large) we will use here only the first zero of the 
Riemann Zeta function up to 10 decimals ((3 small). 

The upper bounds for the amplitude of the waves above, will be compared with the 
results of the numerical experiments performed for the various cases using (3). 



IV. NUMERICAL EXPERIMENTS 



We now present the results of our numerical experiments which was carried out in more 
cases using the Mobius function in (3) up to n = 10 6 . We calculated until k = 10 6 
or k = 10 9 with a scaling factor of 2500 for the fc-axis. These will be compared with the 
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TABLE II: The "amplitude" of ij) for different values of a and (3 



a (3 The function The amplitude 



2 


2 


. 3 

\kic k \ 


= 1^; 2, 2, |)| 


0.000078 


1 


2 


\k*c k \ 


= |V(fc;i,2,|)| 


0.0000292558 


2 


6 


\k*c k \ 


= h/#;2,6,i)| 


0.0210433 


7 
2 


4 


17 2 1 
|fc4 Cfc | 


= |^;l,4,i)| 


0.008411 


3 


3 




= |V(M,3,|)| 


0.0021562 


4 


4 




= 1^(^4,4,1)1 


0.00984936 


2 


4 


II 2 1 

|fc8c fc | 


= l^(fc;2,4,i)| 


0.0052445 



upper bound for the amplitude of the waves of Section 3. The general situation is that for 
moderately values of k (until some tausend) the wave given by the experimental results start 
with zero amplitude, after a minimum with a negative value, increases and seems to stabilize 
at large values of k with oscillations displaced at larger and larger wavelength (proportional 
to log(fc)) and with an amplitude which seems to saturate to a constant value (given in a 
good approximation) by the upper bound (19). Below (Figure 1) we first give the plots of 
the wave for the Riesz case (a = (3 = 2). As remarked in jj], the first intensive calculations 
with very high precision up to k = lOO'OOO (by K. Maslanka) and up to k = 200'000 (by M. 
Wolf) indicated the appearance of oscillations with the first one in the region k = 20'000. 
Our results obtained with (3) confirm for such values the asymptotic limit for the wave with 
an amplitude in agreement with the bound obtained above (A = 0.000078). Notice that the 
minimun for low k values is 0.4 in absolute value as found in a previous work, is much bigger 
then A, which concerns only the asymptotic region of the wave, thus no disagreement! 

Figure 2 and Figure 3 concern two cases of special interest since the decays are expected 
to be the same as for the Hardy-Littlewood case and for the Riesz case. In both cases 
there is agreement with the bound A = 0.0210433 and A = 0.008411 given above but the 
amplitudes are respectively 1000 and 100 time bigger than in the former cases. 

5 7 

In Figure 4 and Figure 5 we give the plots of /cecfc(3,3) and k&c k {AA) where the ampli- 
tudes are found to be in agreement with the theoretical upper bounds given above in Table 
2, too. 

The next special case is the one with a = 2 and (3 — 4. Again, the experimentally 
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3 

FIG. 1: The wave k*Ck for the Riesz case a = f3 = 2 
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FIG. 2: The wave k^cu for the case a = 2,(3 = 6 
detected amplitude agrees well with the theoretical bound given above, i.e A = 0.0052445. 

3 

As a further illustration we compare the wave k&Ck with the asymptotic approximation 
wave given by the wave ip(k; 2,4, |) (case a — 2,j3 — 4). In the range for k from 1 • 10 6 to 
10 • 10 6 the two waves appear to be walking close together arm in arm (Figure 7). Notice 
that in the approximation we considered only the contribution of the first zero given by 
z = \ + HA. 134725141 which appears dominant for low values of (5. 

V. THE CASE a = \ AND (3 INCREASING 

For a = | we will now present the plots of the waves for an increasing sequence of (3 
values i.e. 4, 8, 12, and 20 (in order to investigate the "infinite beta limit" already introduced 
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FIG. 3: The wave k^Ck for the case a = ^,(3 = 4 
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FIG. 4: The wave keck for the case a = (3 = 3 
in our previous work |i7]). We will compute the function 

^~,(3,~) = kh k ( 7 ~,f3) (20) 

which will also be compared with the expression given by the Baez-Duarte formula (19) in 
the asymptotic region k — >• oo. Here we will take into account only the contribution of the 
groundstate of the spectrum i.e z — \ + zl4. 134725141. It is then convenient to introduce 
the new variable x = \og(k). This allow us to control more efficiently the wavelength and 
the amplitude of the wave in the region to be considered (x runs from 8 to 22, so k up to 
3.6 ■ 10 9 ). 

In the Figures 8-11 we present our numerical results for increasing (3 values, which we 
call the "strong coupling limit". 

At the same time it is seen that in this case \ck\ itself is smaller than (c*. is not the critical 
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FIG. 5: The wave k»Ck for the case a = (3 = 4 




-o.oi 

-0.015 



3 

FIG. 6: The wave k&Ck for the case a = 2,(3 = 4 

function!): 

— - 1 S 0.11247897 (21) 

c(i) a 

at least for the case (3 = 4 as already discussed in our previous work [2| concerning only very 
low values of k. Figure 12 confirm this behaviour also for large value of k. For this example 
the region of annihilation of the "eincoming" wave extends up to larger and larger values of 
k. It should be noted that for the critical function i/j (20) the situation is more delicate since 
the value of a possible bound on ip depends on [3 as it is been from our numerical results. 

VI. CONCLUSIONS 

In this work we have extended the numerical experiments obtained in our previous work 
and considered many cases of "waves" (with two parameters a and (3). The experiments allow 
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FIG. 7: The wave k&Ck (lowest curve) and the approximation tp (highest curve) 
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FIG. 8: The wave fc^c^ (lowest curve) and the approximation tp (highest curve), (3 = 4 

to obtain up to 6 oscillations at low (3 values for (20) whose amplitudes are in agreement 
with those given by an extension of an asymptotic formula due to Baez-Duarte for the case 
a = (3 = 2. 

In the process of increasing the values of (3 (at least in the case a — |) it has been 
shown that the width of the "annihilation" region increase with (3 where the amplitude of 
the wave seems (still) to remain bounded. In this connection we may argue that the results 
of our numerical experiments indicate that RH may barely be true due to the behaviour of 
ip in the large (3 limit. It is also conjectured that an absolute bound on ■0, for all a and (3 
allowed, should be given by l^rq^y — 1| — 1-68. In fact Littlewood has shown that on RH, 

^2 J^L is convergent for all e > 0. In the light of the results of the experiments obtained 

n=l nI ' 

so far and with the understanding that our remark is speculative, we belive that even with 
more sofisticated experiments it will be very hard to obtain values of the critical function 
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FIG. 9: The wave k»Ck (lowest curve) and the approximation ip (highest curve), (3 = 8 
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3 

FIG. 10: The wave k^Ck (lowest curve) and the approximation ip (highest curve), (3 = 12 

k a Ck which in absolute value will be greater at large k then those of the infinite (3 limit, as 
commented in the Appendix. 




APPENDIX A 



In the context of the numerical experiments performed so far, it is helpful to obtain a 
crude inequality concerning a bound on the critical function. This is simply obtained by 
setting = 1 in (3). The critical function in the representation of ^ in terms of the 

two parameter Pochammer's polynomials is given by: 

MM 

n=l ^ ' 

where ./V is the maximum value of the argument in the Mobius function considered in a ideal 
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FIG. 11: The wave k^ock (lowest curve) and the approximation ip (highest curve), (3 = 20 




FIG. 12: c k (l,P) for = 4,8, 12, 16,24 (from left to right) 

numerical experiment (TV finite). We then have, introducing the variable x = log (A;) that: 

\h(a,P,N)\ < e^(C(a) - l)e log(1 -^ 

For large N we have: 

\f k (a,(3,N)\<(C(a)-l)e^ x -^ eX 

As an example we consider our case a — | and (3 — 4. Remembering that from Table 
2 the amplitude calculated only with the first non trivial zero is about 0.008411, we may 
ask: for what N and k, \fk(ct, (3, N)\ is bounded by the value 0.008411? For example the 
inequality is satisfied for the followig pairs: 

N = 1000 x > 31 
iV = 10 6 x > 60 
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As a second example we consider the Riesz case (a — f3 — 2). From the Table 2, the 
amplitude (still restricting to the contribution of the first zero) is 0.000078. The inequality 
is satisfied as follows: 

N = 1000 x > 17 
N = 10 6 x > 31 
iV = 10 9 x > 87.2 



Now, still for the Riesz case 



C(*) 



< ^ 1^(5,2,2)^(2,2)1 



fc=0 



assuming as seen in the experiments that k 2 f J Ck(2, f3), i.e. the critical function is bounded 
then we would have that: 

1 



C(«) 



k=0 



Thus, would be different from zero for JH(s) > | + 5, 5 > 0. It should be said that for 
a — | and /3 = 1, in the numerical experiments, the wave seems to decay at zero after a 
few of oscillations, a situation very different from the case a — \ and (5 — 4 (but the above 
inequality apply also in this case). 

So, if the critical function for any (3 is in absolute values bounded for k > K by the value 
of the infinite (3 limit 1.68. . . , the RH should be true. 

In the plot of ip for the case ol — \ and (3 — 4 (Figure 13) up to k = 10 billions we see 8 
oscillations. 
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-0.4 




FIG. 13: k*c k for a = \ and (3 = 4 until k = 10 • 10 9 

Final comments. 
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1. Further application of the crude inequality considered in this Appendix indicates 
that in a ideal experiment using (3), with N = 10 10 , the amplitude of the wave 
for log (A;) < 23 will change at most 10~ 6 time the value 0.008411... obtained with 
N = 10 6 in (3) and compatible with the Baez-Duarte amplitude using only the first 
nontrivial zero of Zeta. This indicates some stability of the numerical experiments in 
the intermediate range log(fc) < 23 (see Figure 13). 

2. In the general case of an ideal experiment with large N in (3) the crude inequality 
indicates also that for k > K(N), infinite, the amplitude of the critical function ip is 
bounded in absolute value by — 1| — 1.68, choosing TV as needed. In the same way 
we may argue that the amplitude of the critical function may be obtained as small as 
we want (in particular smaller then 0.008411...) choosing N as needed for values of 
k > K(N), A - still finite. 

3. One of the open questions is now the following: the critical function at large value of k 
is growing, stabilizing to a "periodic pure wave" with constant amplitude or decaying 
with a zero amplitude? From the results of our numerical treatment we are more in 
favour of the last two cases. 
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